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ABSTRACT 

The  winter-time  transport  of  geostrophically  computed 
angular  momentum  at  500  mb  for  the  period  1  December  1955 
to  28  February  1956  is  obtained  for  Sault  Ste.  Marie, 
Michigan.   This  transport  is  subjected  to  a  cospectral 
analysis  to  compute  the  contribution  of  harmonic  waves  with 
periods  1-30  days  to  the  total  angular  momentum  transport. 
Computations  are  made  for  two  other  stations,  one  east,  and 
one  south  of  Sault  Ste.  Marie  for  comparison  of  transport  as 
related  to  general  circulation  processes. 

The  author  is  deeply  indebted  to  Professor  F.  L.  Martin, 
Department  of  Meteorology,  U.S.  Naval  Postgraduate  School, 
for  his  encouragement  and  invaluable  aid  in  the  preparation 
of  this  manuscript. 
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1.   Introduction 

Vidger  [lOJ  in  his  study  of  the  flow  of  absolute 

angular  momentum  in  the  atmosphere  stated: 

Although  even  a  complete  knowledge  of  the  angular 
momentum  transfer  in  the  atmosphere  cannot  by  itself 
furnish  a  solution  of  the  problem  of  the  general 
circulation,  any  proposed  scheme  for  this  circulation 
should  include  a  means  for  securing  the  angular 
momentum  flows  which  are  observed.   For  this  reason, 
it  is  essential  that  all  possible  observational  knowl- 
edge concerning  these  processes  be  obtained. 

Many  authors  in  recent  years,  Vidger, ' Starr,  and 
Mintz  [10,  7,  3]  for  example,  have  provided  studies  in 
connection  with  the  relationship  of  the  absolute  angular 
momentum  transport  to  the  general  circulation,  or  the  global 
balance  of  angular  momentum.   In  general  these  studies 
provide  quantitative  findings  of  the  spatial  distribution  of 
the  total  momentum  transport  and  establish  the  fact  that 
there  is  a  mean  poleward  transport. 

This  paper,  in  hopes  of  contributing  to  the  knowledge 
of  angular  momentum  transport,  has  as  its  purpose  a  cospect- 
tral  analysis  of  the  angular  momentum  transport  for  waves  of 
periods  1-30  days.   Due  to  the  limited  time  available  for 
the  study,  it  was  confined  to  the  500-mb  level  at  a  station 
in  the  mid-latitudes.   The  primary  station  employed  was 
Sault  Ste.  Marie,  Michigan. 

Several  other  cospectra  are  included  for  other 
stations  for  the  purpose  of  comparison.   However,  the 
spectrum  for  these  comparative  stations  are  for  periods 
2-30  days. 


2.   Mathematical  Formulations 

Widger  [lOJ  has  shown  that  the  rate  at  which  the  ab- 
solute angular  momentum  of  the  zonal  wind  is  increasing  in 
a  portion  of  the  atmosphere  bounded  by  the  earth's  surface, 
the  top  of  the  atmosphere  and  extending  from  the  pole  to  a 
given  latitude  <f>   may  be  expressed  as 

where  dt/  is  any  element  of  volume  within  the  polar  cap  Just 
described,  J5  an  element  of  the  surface  of  the  volume,  ol<3* 
the  projection  of  an  element  of  the  earth's  surface  area  on 
the  meridional  plane  passing  through  a  given  point,  p    is 
the  pressure,  r  the  radius  of  the  earth,  and  "V*  the  east- 
ward frictional  stress  acts  only  on  the  air  at  the  earth's 
surface. 

In  the  first  term  on  the  right  \Jn   is  the  meridional 
component  tr  of  the  velocity  into  the  polar  cap.   Since  the 
vertical  velocities  at  the  earth's  surface  and  top  of  the 
atmosphere  may  both  be  considered  zero,  this  quantity  as  a 
whole  represents  the  rate  at  which  angular  momentum  is 
brought  into  the  region  by  air  motion  across  the  boundary. 
The  second  and  third  terms  on  the  right  represent,  respect- 
ively, the  rate  at  which  angular  momentum  is  changed  due  to 
the  torque  exerted  by  pressure  differentials  across  mountain 
ranges  and  the  westward  torque,  TT^.  ,  exerted  by  frictional 
drag  at  the  earth's  surface.   However,  since  this  study 
will  consider  only  the  transport  through  a  unit  area  of  the 


latitude  wall  at  the  500-mb  surface,  the  terms  involving 
pressure  and  friction  torques  may  be  neglected.  No  attempt 
is  made  herein  to  determine  the  left  side  of  equation  (1) 
from  the  computed  transport  at  500  mb. 

The  absolute  angular  momentum  of  a  gram  of  air  at 
latitude  p    may  be  written 

/1  =  It  r  CdS   <f)   *     O    >  '"C--^}  (2) 

Furthermore,  since  V„~aT  and  the  element  of  surface  may  be 
taken  as  an  element  of  the  latitude  wall,d5-J^Jz  »  taking 
the  first  term  on  the.,  right  in  equation  (1),  the  transport 
of  absolute  angular  momentum  may  be  written  as  the  sum  of 
two  terms, 

rcos  j>  fj  fa^r  d*  dz  (5) 

equal  to  the  rate  of  advection  of  "relative"  angular 
momentum,  and  the  -ft-- momentum  defined  as 

iir'cos^  fff  as    J-k  J  Z  W 

which  is  equal  to  the  rate  of  advection  of  angular  momentum 
due  to  the  earth's  rotation.   Since  JV momentum  does  not 
lend  itself  to  a  cospectral  analysis  and  will  not  be  com- 
puted in  this  study,  the  chief  emphasis  is  placed  upon  the 
average  angular  momentum  transport  through  a  unit  area  and 
per  unit  time  given  by  equation  (3)«   Therefore,  after 
applying  the  geostrophic  approximation  for  U  and  rf  in 
equation  (3)  and  dividing  by  W,  one  obtains  the  average 
angular  momentum  transport 


where  the  hydrostatic  equation  has  also  been  used  in  (3). 

At  first  glance  equation  (5)  appears  to  be  a  summation 
of  N  terms  around  the  latitude  circle  at  a  distance  A  X 
apart,  and  thus  a  representation  of  equation  (3).   However, 
a  slightly  different  interpretation  was  placed  upon  equa- 
tipn  (5),  in  that  AX  represents  a  grid  distance  of  300 
nautical  miles  centered  on  the  station  for  which  the 
computation  is  being  made,  and  A  ^  and  A 2^  are  the  contour- 
height  increments  defining  the  mean  geostrophic  -wind  com- 
ponents 

in  the  neighborhood  of  the  station.   In  addition, A/  was 
taken  equal  to  AX  in  equation  (6).   Thus  the  interpretation 
of  equation  (5)  actually  employed  here  was  to  regard  it  as 
a  time-average,  over  N   equal  time  intervals,  of  values  of 
relative  angular  momentum  transport  centered  at  the  station. 
The  vertical  increment  for  convenience  was  taken  as  Ap  equal 
to  one  millibar  and  A/^  on  the  right  side  of  equation  (5) 
was  scaled  down  to  one  centimeter  instead  of  the  equivalent 
length  of  300  nautical  miles  used  for  the  grid.   The  units 
of  angular  momentum  transport  in  this  study  will  then  be 

cm3-m6  . 

Equation  (5)  shows  that  the  geostrophic  relative  mom- 
entum transport  depends  upon  the  cross  covariance  between 


Uj  and  A/a.   It  therefore  represents  the  average  geostrophic 
relative  angular  momentum  transported  across  a  section  of 
the  latitude  wall  with  grid  length  &/k~  I  &rr>    and  vertical 
thickness  A^r  /tnh' 

The  mean  crosc  covariance  S{t)   of  two  stationary, 
Gaussian  time-series  U(t)  and  /Wt),  at  lagT,  may  "be  ex- 
pressed by  the  convolution  integral 

S(r)-   pi«  jr  f\(t)«r(t+r)dt  (7) 

o 

Frequently  in  meteorological  applications,  O(T)  is 

"evenized"  by  forming--  the  quantity 

S(r)+  5(-t)  (3) 

2 

Lch  is  independent  of  the  sign  of  lag T .  Henceforth,  it 
will  be  understood  that  the  mean  cross  covariance  has  been 
evenized. 

The  cross-spectrum  (or  simply  cospectruia)  Lh  of  the 
two  scries  l((t)  and  sirft)  is  defined  as  that  part  of  the  mean 
cross  covariance  associated  with  a  unit  wave  band  centered 
at  wave  number  h  in  both  series.   In  symbols  C^  may  then 
be  defined  as 

^m  '   r-»ao       p  (9) 

where  P    is  the  averaging  period  and  9<^*°)  an<^  $"■(*}    sire 
the  Fourier  transforms  of  u(t)  and  ^f(t)   respectively,  de- 
fined as  in 


9hW~T-*«  J  L^  e     at     w--^r1 


(10) 


The  circumflex  indicates  the  complex  conjugate. 

If  one  now  employs  the  definition  of  cross  covari- 
ance  ^(X)   for  an  evenized  function,  as  contained  in  equation 
(7),  the  Faltung  theorem  leads  to  the  Fourier  integral 
expression 

so  that  S(v)  and  C(r>)   are  spoken  of  as  Fourier-transform 
pairs.   Equation  (11)  may  be  inverted  to  solve  for  C*  in 
terms  of  S( r).   Also,  the  resulting  inversion  can  be  spec- 
ialized to  the  case  of  a  finite,  but  long,  averaging  period 
P  ,  giving 

Cm '-  l  J  r  S(r)  ,  o  s  ^  v  d  r      »=*??„  (12) 

A  proof  of  equation  (12)  is  given,  from  first  principles, 
in  Appendix  I . 

When  discrete  data  are  employed  at  time  intervals 
of  At  ,  using  a  maximum  lag  of  *><&*,  the  maximum  period, 
or  fundamental,  that  can  be  resolved  is 

■£,•*  ~  I  mat  (13) 

In  equation  (12),  since  T  appears  as  a  summand,  T  will  be 
written  as 

r  =  J>  At       <fc*«,  Ijl, >  m  (14) 

with  &  taking  on  each  of  the  discrete  values  shown  in 
(14).  oince  the  averaging  period  of  equations  (12)  is 
Prm^t,  equation  (12)  leads  to 


C^j^fV^Sjff^^^]     (15) 

I  ir  h  <£  M  -I 

In  the  case  of  do     and  C^  ,  the  resulting  coefficients  are 
divided  by  two,  since  only  one-half  bandwidth  is  available 
at  n-0  and  n~n*>  .   Combining  equations  (5)»  (7),  and.  (8), 
the  cross  covariance  for  A  lags  is 

sinoe  the  averaging  period  of  equation  (7)  is  now  P '(N-k)^t  , 
The  factor  two  in  the  denominator  of  equation  (16)  results 
from  evenizing. 

The  cospectral  estimates  offered  by  equation  (15) 
are  termed  "raw  estimates".   Blackman  and  Tu:ey  [_l\   have 
suggested  several  smoothing  procedures,  the  simplest  of 
which  involves  application  of  a  moving-weighted  average  def- 
ined by 

C*  =.7SCn.\+,SC»  +  \2SC»+[  (17) 

In  the  case  of  rn^hn   there  is  no  higher  value  to  include  in 
the  average  and  the  weighted  average  is  modified  so  that 


3.   General  Procedures 

The  initial  step  required  establishing  what  data 
were  to  be  used  in  making  computations.   Recalling  that  the 
computed  periods  are  determined  by  the  relation  2~rr\&t/f)     , 
where  ho  is  the  maximum  number  of  lags  at  discrete  time 
intervals  At  ,  it  is  clear  that  a  half -day  lag  interval  re- 
quires 30  lags  to  obtain  a  spectrum  for  periods  1-30  days. 
On  the  other  hand,  a  one-day  sample  and  lag  interval  re- 
quires only  15  lags  to  obtain  a  spectrum  for  periods  2-30 
days.   However,  Blackmail  and  Tukey[lJ  state  that  the  best 
results  are  obtained  when  the  longest  lag  is  about  5-10 
per  cent  of  the  length  of  the  record.   This  would  require 
that  the  record,  if  possible,  consist  of  300  -  600  observ- 
ations when  using  the  500-mb  charts  produced  at  0300  GMT 
and  1500  GMT  daily,  or  half  this  number  when  using  only  one 
series  of  charts.   In  order  to  obtain  this  number  of 
observations,  it  would  require  a  minimum  of  5  months  data 
and  the  spectrum  would  be  influenced  by  seasonal  regimes. 
In  order  to  reduce  the  data-processing  load  and  inter-annual 
influences  as  much  as  possible,  it  was  decided  to  use  the 
observations  from  the  three  winter  months  of  December  1955* 
January  1956,  and  February  1956.   Thus  the  longest  lag  was 
only  17  per  cent  of  the  length  of  the  record  for  this 
particular  season.   However,  no  claim  is  made  that  the 
cospectra  obtained  are  necessarily  representative  of  those 
for  other  winter  seasons. 

Lata  were  taken  from  two  sources,  depending  on  the 
station  being  considered.   Since  heights  of  both  the  0300 
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GMT  and  1500  GMT  500-mb  level  were  required  at  Sault  Ste. 
Marie,  Michigan,  the  primary  station  studied  in  this  paper 
facsimile  charts  transmitted  "by  the  U.  S.  Weather  Bureau 
were  utilized.   In  the  case  of  the  other  stations  used  for 
comparison,  the  smallest  period  computed  was  only  two  days 
and  the  data  used  were  the  1500  GMT  charts  published  in 
Daily  Series  Synoptic  Weather  Maps  [8  J. 

From  the  above-mentioned  charts,  the  geopotential 
height  differences  were  obtained  with  a  5*0  degree  (lati- 
tude) grid.   Assuming  that  normal  synoptic  smoothing  was 
applied  in  preparation  of  the  charts,  a  grid  of  this  dimen- 
sion should  not  filter  out  significant  short  waves  that 
have  not  already  been  removed  by  synoptic  smoothing. 

In  five  cases,  where  the  0300  GMT  facsimile  charts 
were  missing,  geopotential  height  differences  were  computed 
from  the  actual  observed  wind  listed  in  Northern  Hemisphere 
Data  Tabulations  L9J«   The  observed  winds  were  reduced  to 
geostrophic  values  by  subjectively  compensating  for  the 
degree  of  cyclonic  or  anticyclonic  curvature  associated  with 
the  flow. 

The  smoothed  cospectra  of  geostrophic  relative 
angular  momentum  transport  were  obtained  by  using  these 
data  in  equations  (15),  (16),  (17),  and  (18),  and  are  graph- 
ically displayed  in  figs.  3  through  6.   All  computations 
except  for  the  final  smoothing,  utilizing  equations  (17),  (18), 
were  accomplished  on  the  CDC  -  1604  digital  computer* 


4.   Limitations  and  Sources  of  Error 

The  computational  results  of  this  study  have  been  based 
on  discrete  data  using  equal  intervals  A"tr  .   In  one  series 
of  computations,  indicated  by  the  solid  line  in  fig.  3,  the 
data  interval  was  reduced  to  12  hours  in  order  to  determine 
the  effects  of  sampling  fluctuations  and  of  atmospheric 
oscillations.   In  all  other  cospectral  analyses,  the  data 
interval  was  24  hours.   Throughout  the  study  the  maximum  lag 
was  15  days. 

The  period  associated  with  the  smoothed  cospectral 
estimate  is  Z\m&*-/r\   with  Kl^i,  2,  •*'  '  t  ho  .   On  the  other  hand 
the  central  frequency  £,  of  the  frequency  band  is 

f»  ~     2 * *t  hr''^ '* 

and  the  highest  frequency  resolvable  occurs  for  h  =  >d  ;  this 

frequency  is  called  the  Nyquist  frequency  T^  .   The 
cospectral  estimate  corresponds  to  the  frequency  band 

(»-  Yt)/zn/$t    to     (h-f  [/z)Jt 

It  should  be  noted  that  it  is  difficult  to  resolve 
frequencies  contributing  to  the  cospectrum  at  the  longer 
periods.   For  example  when  h  -  t   ,  the  frequency  band  spans 
the  range  */ij   cycle  /  day  to  l/fco   cycle  /  day. 

In  order  to  assess  the  reality  of  sharp  peaks  or 
troughs,  Panofsky's  criterion  f or  X/V  distribution  of  the 
power  spectrum  has  been  followed.   Strictly  speaking  a  much 
more  rigorous  method  of  assessing  significance  would  have 
been  achieved  by  testing  the  cospectral  correlation 
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function  defined  by 


R    .  - £sa_ 


'J' 


Here  C>>>  is  the  cospectral  estimate  for  the  in   *    harmonic, 
and  Pk),^  and  r^^are  the  power  spectral  estimates  for  U. 
and  at  at  the  same  harmonic.   However,  the  power  spectra 
computations  were  not  available  in  time  for  the  completion  of 
this  aspect  of  the  work.   Hence,  the  more  qualitative  ap- 
proach indicated  in  fig.  1  was  attempted  with  some  misgiv- 
ings.  It  must  be  emphasized  that  when  the  term  "95^ 
confidence  limit"  is  used,  it  is  only  by  analogy  with  corres- 
ponding phraseology  to  be  found  in  Panofsky  and  Brier  |_4  , 
p.  14-5  J  •   The  actual  confidence  may  be  at  a  somewhat  lower 
level. 

Applying  the  X/y   distribution  to  the  co spectrum,  a 
peak  or  trough  covld  not  have  happened  randomly  (at  the  95% 
level  of  belief)  provided  the  peak  or  trough  in  the  cospectral 
estimate  lies  outside  the  range  £>P  ,  defined  as 

aPhpXxVv)  (19) 

where  JT    is  the  mean  cospectral  intensity  between  the  peak 
and  trough,  X"QS/y   is  the  chi-square  value  at  the  .05  level 
for  V       degrees  of  freedom  (^  V  -  c.  Nf  w\  "  /i,  iz  12.  J   .  Accord- 
ing to  Panofsky  and  Brier  ^4-,  p.  55J  the  limiting  value  of 
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Fig.  1.    Illustration  of  the  X/\>    test  of 

significance  at  the  95%  level  of  belief. 


A  second  source  of  error  which  may  considerably  bias 
the  accuracy  of  cospectral  estimates  is  a  phenomenon  called 
"aliasing".   This  phenomenon  removes  energy  from  higher 
frequencies  than  those  sampled,  and  may  convert  it  to  lower 
frequencies.   Fig.  2  illustrates  how  aliasing  may  occur. 


Fig.  2.    Illustration  of  the  effects  of  sampling 
variations  in  giving  an  aliased  spectrum, 
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Suppose  a  periodic  function,  such  as  that  illustrated 
in  fig.  2,  exists  for  "both  U^(t)  and  W(t)   ,  and  consider  the 
period  of  the  oscillation  shown  to  be  2At  „   If  lags  are 
taken  at  intervals  of  3  At:  ,  one  is  investigating  the 
cospectral  estimate  at  a  period  of  Q/^'t  .   if  the  sampling 
times  corresponded  to  points  A,  B,  and  C,  one  obtains  a 
rather  large  "aliased"  cospectral  estimate  for  the  period  Q At 
However,  if  the  sampling  points  had  been  E,  F,  and  G,  the 
aliased  cospectral  contribution  at  period  G^t  would  be 
zero.   Note  that  the  oscillation  could  be  a  real  atmos- 
pheric oscillation  or, meteorological  "noise"  of  an  obs- 
ervational character.   The  primary  atmospheric  oscillations 
of  periods  shorter  than  one  day  are  the  semi-diurnal  (12-hr) 
wave  and  a  10.5-hr  wave,  originally  described  by  Pekeris  £5] 
However,  the  amplitude  of  these  waves  is  quite  small. 

Blackman  and  Tukey  (.1J  have  shown  that  the  frequencies 
which  are  the  principal  aliases  of  frequency  f„  are  Z'f^'tf^ 
^i^f*  \  Qf^tf    i  etc.,  where  ^-'/^tis  the  Nyquist  freq- 
uency.  However,  the  effect  of  aliasing  is  most  pronounced 
at  frequencies  T^    near  the  Nyquist  frequency.   Griffith  et  al 
\Z J  suggest  that'  one  way  of  eliminating  aliasing  from  the 
cospectrum  is  to  delete  the  cospectral  estimates  for  the 
last  20  to  4-0  per  cent  of  the  higher  frequencies.   However, 
an  experiment  was  performed,  using  data  for  Sault  3te.  Marie, 
to  assess  the  extent  of  aliasing. 

Fig.  3  shows  tliree  cospectra  for  Sault  Ste.  Marie 
based  on: 

(a)   the  1500  GMT  reports  at  24~hr  intervals  (curve  A), 
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(b)  the  0300  GMT  reports  at  24-hr  intervals  (curve  B) , 
and 

(c)  "both  the  0300  GMT  and  1500  GUT  reports  at  12-hr 
intervals  (curve  G). 

It  may  be  noted  from  fig.  3  that  all  cospectral  est- 
imates at  periods  greater  than  2.73  days  are  in  fairly  close 
agreement.   However,  curves  A  and  B  deviate  rather  markedly 
at  their  Nyquist  frequency  of  Vi   cycle  /  day. 

In  the  case  of  curve  C  there  should  be  little  or  no 
aliasing  at  periods  greater  than  two  days  according  to 
Griffith's  criterion,  since  50  per  cent  of  the  spectrum 
lies  between  the  periods  one  and  two  days.   On  the  other 
hand,  at  all  other  stations  it  will  be  assumed,  henceforth, 
that  all  cospectral  estimates  at  periods  2.73  days  or  larger 
are  relatively  unaffected  by  aliasing. 
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5*   Interpretation  of  Cospectra 

Figs.  3  through  6  present  the  six  cospectra  computed 
for  this  study.   In  these  figures  the  ordinate  is  the  co- 
spectrum  and  the  abscissa,  plotted  for  convenience  on  a 
logarithmic  scale,  is  the  period. 

In  assessing  the  reality  of  the  various  peaks  observed 
in  fig.  3>  the  procedure  embodied  in  equation  (19)  will  be 
applied  only  to  curve  C,  as  the  other  two  curves  do  not 
deviate  significantly  except  at  a  period  of  3.75  clays,  which 
will  be  discussed  separately.   In  order  to  compute  a  mean 
cospectral  estimate  for  the  peak  at  a  period  of  10  days,  it 
was  assumed  that  the  cospectral  estimate  at  30  days  is  the 
apex  of  a  trough.   Inspection  of  the  95%  and  5%  fiducial 
limits  established  from  the  mean  cospectral  values  on  each 
side  of  the  peak  and  trough,  indicate  that  the  peaks  and 
trough  associated  with  periods  10  and  6  days,  respectively, 
are  significant.   The  peaks  and  troughs  have  only  a  5% 
probability  of  randomly  exceeding  these  limits;  therefore 
both  are  considered  to  be  real. 

In  the  case  of  curve  A,  the  peak  appears  not  to  be 
significantly  different  when  compared  to  the  fiducial  limits 
computed  for  curve  G.   However,  limits  computed  from  mean 
values  of  curve  A,  although  not  shown  here  in  order  to  avoid 
confusion  with  those  of  curve  G,  show  the  value  at  the  peak 
on  curve  B  is  considerably  greater  than  the  limit  to  the 
left  side  and  slightly  greater  than  the  limit  to  the  right 
side,   although  the  peaks  at  3.75  days  are  not  markedly  sig- 
nificant from  the  limit  on 'the  right  side,  there  is 
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justification  for  stating  these  peaks  seem  to  be  real  when 
one  considers  both  the  left  and  right  side  limits.   Using 
the  significance  test  discussed  above,  the  reality  of  any- 
minor  peaks  in  fig.  3  at  periods  of  3  days  or  less  is 
doubtful. 

The  curve  in  fig.  4-  is  the  average  cospectrum  for 
Caribou  during  the  period  1  December  1955  to  28  February 
1956.   A  one-day  interval  between  data  points  was  used  to 
compute  this  cospectrum.   An  inspection  of  the  fiducial 
limits  shown  in  fig.  4-  indicates  that  the  peak  at  6  days  is 
significantly  different  from  the  limits  on  both  the  left  and 
right  side  of  the  peak,  and  it  is  therefore  considered  to  be 
real.   As  to  the  small  peak  at  approximately  3.3  days,  the 
confidence  interval  shown  to  the  left  is  quite  representative 
for  this  peak.   It  can  therefore  be  concluded  that  this  peak 
(at  3.3  days)  is  probably  not  real,  as  it  does  not  lie  beyond 
the  confidence  range.   However,  it  can  not  be  inferred  from 
this  that  the  peak  is  a  result  of  sampling  variations,  since 
the  test  does  not  with  certainty  disprove  the  reality  of  the 
peak. 

The  curve  in  fig.  5  is  also  for  Caribou,  but  was  com- 
puted from  data  taken  during  the  winter  previous  to  that 
used  in  fig.  4.   The  limits  in  this  case  show  the  trough  and 
peak  at  7.5  and  5  days,  respectively,  to  be  real.   The  trough 
at  3.75  days  is  questionable,  since  it  is  not  significantly 
different  from  the  limit  to  the  right. 

One  of  the  more  striking  features  in  each  of  the 
cospectra  is  the  large  transport  associated  with  the  longest 
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period  and  the  fact  this  transport  has  the  same  sign  as  the 
total  transport. 

In  general  it  is  difficult  to  associate  a  particular 
period  of  the  cospectrum  with  a  specific  space  wave,  since 
the  period  is  a  function  of  both  the  wavelength  and  speed 
of  the  space  wave.   It,  thereby,  seems  possible  for  a  trans- 
port at  a  specific  period  to  have  contributions  from  both 
long  and  short  space  waves.  This  would  be  the  case  if  the 
ratio  of  their  wavelengths  is  equal  to  the  ratio  of  their 
speeds.   However,  assuming  that  the  Rossby  wave  equation 
gives  a  valid  relationship  between  the  wave  speed  and  wave- 
length for  any  given  mean  zonal  wind,  it  may  be  expected 
that  long  space  waves  move  eastward  slowly  and  short  ones 
more  rapidly.   In  view  of  this  assumption,  the  transport 
associated  with  a  period  of  thirty  days  may  be  attributed 
to  the  contributions  from  long  space  waves. 

Furthermore,  large  values  of  transport  by  the  longest 
period  were  to  be  expected,  for  an  examination  of  the 
monthly  mean  500-mb  charts  of  December  1955  through  1956 
showed  3ault  Ste.  Marie  and  Caribou  were  on  the  western 
side  of  a  mean  long  wave  trough,  hence  the  net  transport 
would  most  likely  be  southward  (negative).   However,  for  the 
second  cospectrum  at  Caribou  (fig.  5) »  the  period  December 
1954-  through  February  1955  was  selected  by  ascertaining 
through  monthly  mean  500-mb  charts  that  Caribou  was  east  of 
the  trough  and  a  northward  (positive)  transport  could  be 
expected.   In  the  case  of  Athens,  Georgia,  the  main  reason 
for  selecting  this  station  to  provide  a  comparative 
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cospectrum  is  the  fact  it  is  almost  due  south  of  3ault  Ste. 
Marie.   One  could  expect  these  two  stations  to  have  a  net 
transport  of  the  same  sign,  but  examination  of  the  monthly- 
mean  500-mb  charts  show  the  mean  long  wave  trough  to  have  a 
southwest  to  northeast  tilt  of  sufficient  magnitude  to  place 
the  trough  west  of  Athens.  As  a  result  the  net  transport  for 
Athens  was  northward,  while  that  of  Sault  3te.  Marie  was 
southward.   The  indicated  convergence  of  momentum  transports 
should  then  act  in  such  a  way  as  to  increase  the  mean  zonal 
westerlies.   Note  also  that  the  only  pronounced  peak  in  the 
cospectrum  at  Athens  occurs  at  a  period  of  30  days,  and  that 
the  cospectral  values  diminish  with  decreasing  period.   This 
is  probably  due  to  the  fact  that  short  period  waves  do  not 
penetrate  appreciably  into  this  latitude  belt. 

Interpretation  of  significant  transports  associated 
with  shorter  periods  requires  reference  to  the  earlier 
statement,  that  short  waves  are  likely  to  move  rapidly.   In 
view  of  this  concept  the  momentum  transport  evident  near 
periods  of  four  days  is  attributed  to  waves  of  approximate 
space-wave  numbers  6  to  10,  which  Saltzmann  and  Fleisher 
L6J  have  called  "cyclone  waves".   However,  the  momentum 
transport  peaks  at  periods  near  10  days  are  presumably  due 
to  long  waves  of  space-wave  numbers  1  to  5« 
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6.   Conclusions 

It  nay  be  concluded,  tentatively,  from  the  cospectra 
presented  in  figs.  3  through  6  and  the  preceding  discussion 
that: 

(a)  Haximum  transport  of  geostrophic  relative  angular 
momentum  is  accomplished  by  disturbances  with  the  longer 
periods. 

(b)  The  peaks  with  periods  near  10  days  are  related 
to  Ion,  space  waves, 

(c)  The  direction  of  the  net  transport  is  highly  de- 
pendent on  the  position  of  the  mean  long  wave  trough 
relative  to  the  observing  station. 

(d)  Cyclone  waves  of  space-wave  numbers  6  to  10  may 
contribute  to  the  maximum  momentum  transports  associated 
with  periods  near  4-  days. 

(e)  The  maximum  transports  at  periods  of  30  days  may 
be  attributed  to  the  circulation  index  cycle. 
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Fig.  4  The  cospectrum  of  the  mean  relative  angular  momentum 
transport  at  the  500-mb  level  over  Caribou,  Maine  for 
the  period  1  December  1955  through ' 28  February  1956. 
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Fig.  5  The  cospectrum  of  the  mean  relative  angular  momentum 
transport  at  the  500-mb  level  over  Cs  Lbou,  Mad 
for  the  period  1  December  1954-  through  28  February  1955< 
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Fig.  6  The  co-spectrum  of  the  mean  relative  angular  momentum 
transport  at  the  500-mb  level  over  Athens,  Georgia 
for  the  period  1  December  1955  through  28  February  1956. 
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APPENDIX  I 
Derivation  of  the  Cospectral  Function 

The  cross  covariance  at  lag  T  ,  of  two  functions  of 
time,  is  given  by 

S(r)  rrtt  t)0  u  (t) .  W(t  f  r )  d  fc 
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The  first  cosine  term  integrates  to  zero  for  all  cases 
m,n  /  0.   The  second  cosine  term  integrates  to  zero  for 
all  cases  m  /  n.   When  m  =  n,  integration  gives 
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The  evenized  cross  covariance  equals  one-half  the  sum 
of  the  cross  covariance  at  T  and-T  , 

Defining  the  unsmoothed  cospectrum  as 
Ch^^co.-U-  &„) 
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the  Fourier  transform  of  C(n)  is  obtained  as 

o 

Multiplying  both  sides  by  Co% — y= —  a  t    and  integrating  from 

0  to  T  gives 

4TThT 
For  all  n  /  C,CoS = —  integrates  to  zero,  giving 

where  P    is  the  integration  limit  and  T  is  the  fundamental 
period. 
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